How loud can Schwarzschild black holes ring? 
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A numerical procedure is described for the maximization 
of the energy diffusion due to the backscattering of the gravi- 
tational radiation. The obtained maxima are solutions domi- 
nated by low frequency waves. They give rise to robust grav- 
itational ringing, with amplitudes of the order of the original 
signal. 
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I. MOTIVATION 

The detection of gravitational radiation by the first 
generation of laser interferometers requires from theorists 
and numericists developing of the so-called templates for 
the gravitational waves (GW) produced by binary black 
holes (BH) jjj , p| . That can be done for the inspiral of 
BH while there are significant problems as it concerns 
the merger phase of BH , where no numerical data are 
available as yet. In this context the next, ringing phase, 
of the collapse of binary BH is interesting because GW 
produced during this evolution period are independent 
to a degree on the details of its preceding history. Nu- 
merical simulations demonstrate that irrespective of the 
initial data, during some "intermediate" time the outgo- 
ing GW are dominated by a (fundamental) quasinormal 
mode (QM) B, whose oscillation period and the damp- 
ing coefficient depend only on characteristics of the final 
black hole. It is not clear, however, how energetic || are 
the QM and how big is the ringing amplitude. The full 
answer to this question would require detailed informa- 
tion on the former evolution of the system, which is not 
available at present. This paper aims at finding "upper 
limits" , that is identifying the most favourable situations 
for ringing and comparing their strength - amplitudes 
and the energy content - with the initial energy and am- 
plitudes. The amplitudes of strongest ringing modes are 
of the order of the amplitudes of GW that generate them. 



II. EXTREMIZING THE DIFFUSED ENERGY 

A. Prerequisities. The space-time geometry is defined 
by a line element, 
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ds 2 = - 1 )dt 2 -i — 
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dR 2 + R 2 dil 2 , (1) 



where t is a time coordinate, R is the areal radius and 
dtt 2 = d9 2 + sin 2 0d<p 2 is the line element on the unit 
sphere, < <f> < 2w and < 9 < ir. The Newtonian 
gravitational constant G, and c, the velocity of light are 
put equal to 1. 

We will study the propagation of polar modes \& of the 
quadrupole GW in a Schwarzschild background; they are 
ruled by the Zerilli equation B 
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Here r) R = 1 - 2m/ R, V{R) 

and r* = R + 2mln( - 1). 

There exists a positive conserved energy E M; its 
quasilocal contribution coming from a region of a Cauchy 
hypersurface Sj exterior to a sphere of radius R reads 



E(R,t) = 2tt / drp(r,t). 



(3) 



Here p is the reduced energy density, p = ((d t ^) 2 + 
{d r *-%>) 2 + V* 2 ) /r, B , The initial data * and <9 t * are as- 
sumed to be purely outgoing, smooth and to be nonzero 
outside a sphere of a radius a > 2m. Thus p is smooth 
and vanishes inside the sphere of radius a. 

Let rWj) be an outgoing null geodesic that originates 
at (R,t). By T(R 0l t ),(R,t) wn l be understood a segment 
ofF^ ^) ending at (R,t). A straightforward calculation 
shows that the rate of the energy change along TfRm is 
given by 



(d t + d r *)E(R,t) 
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Thus E(R, t) = E(a, 0)-2tt f Q * dt 



where the integral is done along ^( a ,o),(R,t)- Some of the 
energy can diffuse inward due to the scattering off the 
curvature of the spacetime H . In the limiting case case 
when the integration contour coincides with rV a) o), the 
energy limit i?s(a,0) = limj^oo E(R(t), t) is the ana- 
logue of the Bondi mass. We will define the diffused 
energy as 5E a = E(a, 0) — Eb{cl, 0). Our aim is to find 
initial data, that maximize the diffusion parameter 
k = 8E a /E(a,0) for a selection of a's. Let us remark 
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that, from our experience, the computational time is pro- 
portional to the square (a/(2m)) 2 ; therefore the max- 
imization procedure would not be feasible in the case 
a » 2m . Fortunately, it is known from analytic esti- 
mates, that for big a's the diffusion is negligible, since 
k < C{2m/a) 2 (where C is of the order of 10) @, Q. 
Therefore it suffices to focus to the range of relatively 
small distances - a being of the order of the Schwarzschild 
radius 2m - and in that case the numerical methods ap- 
pear to be efficient. 

B. Method. It is advantageous to search a solution 
of Eq. ^ in the form * = * + 5 ffl. Here * is a 
known function that is constructed from initial data while 
5 is the unknown part of the sought solution. In explicit 
terms 



* ee * (r* - t) 



*i(r*-t) , * 2 (r*-i) 
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where #j(r* -t),i = 0,1,2 fulfill the relations d t ^i = 
3^0 and dt^ 2 = \&i — mdt^x- $ solves the Zerilli equa- 
tion in Minkowski space-time (m = 0), in which case it 
represents a purely outgoing radiation. The function 8 
satisfies the equation 

(-i? + ^ = 7i + (K-6|)f*o+5 + ^ + 
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with initial data S = dt5 — at t = 0. 

Eq. (0) and the relations between functions imply 
that initial data W and dt^ are entirely determined by the 
specification of m turn is determined (in a region 

outside the cone defined by T a ) by choosing d t ^o (with 
dt^o\a = 0) in an interval (r*(a), oo) and *o = for R = 
a. We have found that the numerical calculations can 
be significantly facilitated if the following procedure is 
applied. First, expand dt^o in an n-element Chebyshev 
polynomial basis (/j) S; one has 



5 t * (r*,t = 0) = Ef =1 C i / i (r*). 



(7) 



The initial energy becomes in this basis quadratic in Ck ■ 
E(a,R,t = 0) ee 2vr f* drp(r,0) = T^^B^a, R)CiC k . 
Second, let \& f h (r* , t) be a solution of (0) with initial data 
being defined by (r* ) . Obviously * / fc (r* , i) depends 
only on the initial data within the null past cone with 
the apex at (r*(a), t). Then the solution of Eq. ph reads 



*(r*,t)=E2 =1 C,* /t (r*,t). 



(8) 



Let Ri = R(r* = t + r*(a)) and R 2 = R(r* = 2t+r*(a)). 
The energy diffusion through r/ Q .o),(i?!,t) is of the form 



The quadratic form 5E a (t,Ri) is nonnegative; since 
E(a,R,t — 0) is positive, the two forms can be simul- 
taneously diagonalized. The task of comparing the ex- 
tremal ratio of SE a and E(a, R, t = 0) can be reduced to 
the search of the maximal eigenvalue of the generalised 
eigenvalue problem AX = XBX. That has been accom- 
plished with the use of the fast EISPACK [Q package. 
The typical step used in numerics was A = 0.02 and n 
was of the order of one thousand, depending (weakly) on 
the size of the initial support. 

C. Extremal initial data. Below the mass m has been 
normalized to unity. The method described above con- 
sists in comparing local energy expressions - the energy 
diffused through a segment Ti a ,o),(Ri,t) is compared with 
the energy content inside the region a < R < R 2 of the 
initial hypersuface. Fixing a and R 2 , one finds dt^Q 2 
(the upper index is put here in order to stress the local 
character of the procedure) and, through (||), initial val- 
ues of the locally extremizing solution ^ R2 . With the 
increase of R 2 , while keeping a fixed, the function dt^Q 2 
changes. In the limit one obtains the sought extremizing 
solution, * = lim/^^oo fy R2 . In the numerical practice 
the integration region must be finite. 

Fig. 1 shows the initial profile dt'i'o for a — 2.1 and 
for various R 2 . 
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1. Function dt^o in dependence on R 2 . The abscissa 



5E a (R,t) = El k=1 A ik (a, R)CiC k . 
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(R— axis) is in the (decimal) logarithmic scale. Here a — 2.1; 
solid line, long broken line and short broken line correspond 
to r*(R 2 ) = r*(2.1) + 6000, r*(R 2 ) = r*(2.1) + 8000 and 
r*(R 2 ) = r*(2.1) + 10000, respectively. Dotted line depicts 
the case with r*(Ra) = r*(2.1) + 16000. 

The dependence of d t ^o on R 2 suggests that 
limfl^oo dt^Q 2 ~ outside some region of compact sup- 
port. The parameter k was extremized in that class of 
initial data which is characterized by the vanishing of 
d t ^o outside some bounded region. Therefore the cho- 
sen * bears on an asymptotically constant value ]l3| . 
Fig. 2 shows initial profiles of dt^o for a = 2.1, a = 3.1 
and a = 4. 

The obtained extremizing initial data have a finite 
total energy. Fig. 3 shows initial energy profiles for a = 
2.1, a — 3.1 and a = 4. It is interesting to notice how 
efficient is the backscatter - k ranges from c. 90% (a=2.1) 



2 



and c. 16 % (a=3.1) to c. 4% (a=4). Notice also that the 
region with the largest contribution to the initial energy 
widens with the increase of a. 
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FIG. 2. Function ftv&o for various values of a. The ab- 
scissa (R— axis) is in the (decimal) logarithmic scale. Solid 
line, long broken line and short broken line correspond to 
a — 2.1, a = 3.1 and a = 4, respectively. 
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FIG. 3. Initial energy E(R) = E(a, 0) - E(R, 0) (y- axis) 
as a function of r* — r*(a) for a = 2.1 (solid line), a — 3.1 
(long broken line) and a — 4 (short broken line). 

III. EVOLVING EXTREMAL INITIAL DATA 

A. Ringing is loud. Our primary aim is to find initial 
data that give the strongest possible ringing. QM are 
evidently born by a subclass of multiple backscattering 
(see an explanation in Q]). The diffusion energy SE a 
bounds the energies of QM, the tail term and also of 
the radiation falling to a black hole. While we do not 
have analytic estimates of the shares of the particular 
contributing terms in SE a , it is obvious that configura- 
tions with large k have some room for robust oscillations. 
For that reason, instead of obeying the commonly used 
method of trial and error, we study GW defined by the 
extremal initial data. In the numerical calculations 
we use the splitting = & + S (see Sec II. B), since the 
numerics is then more precise. 

Fig. 4 presents the radiation corresponding to the 
a = 2.1, a = 3 and a — 4 initial pulses as seen by an "ob- 
server" situated at r* = 280 + r* (a) . The x = point of 
the abscissa corresponds to the moment of time t = 280. 



This train of initial data that moves with the speed of 
light is seen earlier (t < 280) and it lies to the right from 
x = 0. To the left from x — we have t > 280; in the 
absence of the backscattering there would be no signal at 
all. Notice that in the cases a — 2.1 and a — 3.1 the am- 
plitude of the strongest ringing mode is about one third 
of the amplitude of the original radiation, as represented 
by the function >F. The amplitude of the strongest ring- 
ing mode decreases with the increase of a - that is when 
the initial pulse is moved away from the horizon - but in 
all cases it is of the order of the initial amplitude. 
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FIG. 4. Temporal dependence of In |^| (y— axis), as ob- 
served at ?"o = r*(a) + 280. Solid line, long broken line and 
short broken line correspond to a = 2.1, a = 3.1 and a — 4, 
respectively. 

B. Birth and death of ringing modes. QM boil down 
predominantly in the potential valley 2 < R < R cr w 3.1, 
where the maximal point R cr is defined by djiV\ji cr = 0. 
Their complex frequencies carry information about the 
curvature of the background geometry in this region. 
This is why the characteristic features of the ringing 
modes - the period and the damping coefficient of the 
dominant mode - do not depend on initial data. Let us 
point out, however, that the amplitudes of the QM de- 
pend on the frequency profile and amplitudes of initial 
data. We examined QM in various observation points. 
Fig. 5 (with a = 4) shows clearly that the number of 
nodes in a radiation pulse decreases with the distance. 
At the same time the dominant oscillations survive and 
they have the same character (but their amplitude can 
significantly increase) while the tail part extends. 

Fig. 5 shows that there are many oscillations at 
r* = r*(4) + 1, which gradually die when the obser- 
vation point is moved away - only one node is seen at 
r* = 280 (Fig. 4). This observation that we make is 
probably novel in the literature, but its explanation can 
be standard, within the scenario described in J14j|. 
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FIG. 5. Backscattered radiation (values ln|$| put 
on y— axis), a — 4, as seen by an observer located at 
j-q = 1 + r*(4) (solid line), Tq = 10 + r*(4) (long broken 
line) and Vq = 100 + r*(4) (short broken line). 



IV. CONCLUSIONS 

We have found initial data that correspond to maxi- 
mal values of the diffusion parameter k; we believe that 
these are optimal data for having strong QM. The largest 
amplitudes of QM constitute a fraction (c. 1/3) of the 
largest amplitude within the main pulse of the radiation. 
Our earlier experience Jig ] in numerics suggests that the 
effect of the backscatter depends (fixing the distance a) 
mostly on the frequency, and that if the dominant fre- 
quency is low (as compared to l/m), then the diffusion 
parametr k does not depend strongly on details of the 
profile of initial data. This would imply that in the real 
collapse the ringing can be perhaps less vigorous than 
above, but still of interest, with amplitudes smaller by 
perhaps one order than those of the main pulse of GW. 
Similar conclusions concerning the strength of QM are 
valid also for the axial quadrupole GW and for the dipolc 
electromagnetic waves | fl7| . The main difference between 
the polar and axial perturbations is that polar QM are 
stronger and they seem to be more "generic", in a sense 
that will be explained elsewhere jl7). The same analysis 
can be done for higher GW multipoles. 

Referring to a specific case of collapse of BH, we want 
to stress the following. Our results can be applied to 
the head-on collision of two spinless BH. If one uses the 
close approximation limit , then the outgoing part of 
a radiation that was produced during the earlier phases, 
the merger and the inspiral, should be interpreted as our 
initial data. If these initial data are close to " maximal" , 
in a sense used earlier, then the ringing amplitudes could 
be of the order of GW created during the coalescence of 
BH. 
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